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The production of charmonia is investigated for heavy-ion collisions from SPS to RHIC energies. 
Our approach incorporates two sources of J/^P yield: (i) a direct contribution arising from early 
(hard) parton-parton collisions, subject to subsequent nuclear absorption, quark-gluon plasma and 
hadronic dissociation, and (ii) statistical production at the hadronization transition by coalescence 
of c and c quarks. Within an expanding thermal fireball framework, the model reproduces J/^P 
centrality dependencies observed at the SPS in Pb-Pb and S-U collisions reasonably well. The 
study of the ^''/'I' ratio at SPS points at the importance of the hadronic phase for ^' interactions, 
possibly related to effects of chiral symmetry restoration. Predictions are given for the centrality 
dependence of the Nj/q,/Ncc ratio at full RHIC energy. We also calculate the excitation function of 
this ratio. The latter exhibits a characteristic minimum structure signalling the transition from the 
standard J/'i/ suppression scenario prevailing at SPS to dominantly thermal regeneration at collider 
energies. 



I. INTRODUCTION 

The main objective of ultrarelativistic heavy-ion physics, both experimentally and theoretically, is to create and 
study hot and dense matter of strongly interacting constituents. At sufficiently large temperatures and densities, the 
theory of strong interactions (QCD) predicts a phase transition from the hadronic world into a chirally symmetric 
state of free quarks and gluons, the so-called Quark-Gluon Plasma (QGP). However, the transient nature of this new 
state of matter renders its identification very complex. After more than a decade of experiments, a wealth of exciting 
results has emerged from the SPS program as well as from recent data taken at the the Relativistic Heavy-Ion Collider 
(RHIC) at Brookhavcn. However, as of now, no conclusive evidence of QGP formation has been discerned. It has 
rather become clear that the discovery of the QGP would not emerge from a single signature but from a combination 
of several in a comprehensive scenario. 

Among the probes of QGP formation, the J/'i> meson plays a central role. Due to its large binding energy and 
rather small size, it is expected to be only mildly affected by interactions in a hadronic environment. On the contrary, 
as suggested first by Matsui and Satz ^ in 1986, if a QGP is formed, J/^! formation should be inhibited due to 
color charge screening by free-moving quarks and gluons. From a slightly different perspective, a similar effect can be 
attributed to parton-induced destruction In either picture, the J/^* abundance observed in heavy-ion collisions 
is linked to properties of the initial hot and/or dense phases of the produced matter. However, the identification of 
the plasma effect requires a careful assessment of other mechanisms for J/^E" destruction in the course of a heavy-ion 
reaction. Most notably, absorption on high-energetic nucleons at early times induces a substantial suppression of 
typically 60% in central collisions. Although its dependence has been thoroughly investigated experimentally in both 
proton and light-ion induced reactions with various targets ^, |^, ^ theoretical uncertainties in the application to 
the heavy-projectile case persist [Q. Furthermore, despite the naive expectation given above, the current status in the 
description of inelastic charmonium scattering on "comoving" (secondary) hadrons in the late stages of the collision 
is not satisfactorily under control. This renders the identification of the effect of the quark-gluon plasma formation 
on J/^ abundances more difficult to establish. 

Recently, an additional source of J/^! production in heavy-ion collisions has been suggested. Within thermal 
model frameworks ||, 0, which successfully describe light-hadron ratios, statistical production of charmonium 
at the QCD hadronization transition has been evaluated |lj, |l2|, Gazdzicki and Gorenstein found that 

statistical production alone can account for the observed centrality dependence of J/'i' yields at full SPS energy, 
thereby deducing a hadronization temperature consistent with standard chemical freezeout. However, no reference 
was made to a dynamical origin for cc creation, nor to open-charm states. Braun-Munzinger and Stachel ^ 
started from the dynamically well justified JTst hypothesis that at SPS energies all cc pairs are exclusively created in 
primordial (hard) parton-parton collisions. Open- and hidden-charm hadrons are then populated at the hadronization 
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transition according to their thermal weight including a charm-quark fugacity to match the primordial cc abundance. 
With the latter taken from extrapolations of N-N collisions, a substantial fraction of J/^ mesons is found. In 
particular, the VP' over ^! ratio, which in Pb-Pb reactions has been observed to rapidly approach its thermal value 
for Npart > 100, is well described. Along similar lines, Gorenstein et al. aim at explaining the J/^E* yields at 
SPS solely in terms of statistical production. They conclude that an open-charm enhancement factor of up to ~ 6 in 
central Pb-Pb collisions (w.r.t. the standard value inferred from TV- TV collisions) is required to optimally reproduce 
the NA50 data Finally, Thews et al. assessed cc coalescence in an expanding QGP fireball by solving 

rate equations for J/W+g <-> cc reactions (gluon-induced "photo" -dissociation and its reverse reaction), with the J/^ 
binding energy assumed to be at its vacuum value at all times in the evolution. 

All thermal approaches share the common feature that, as first pointed out in Ref. ||l2|, at RHIC energies a copious 
production of charm-quark pairs {Ncc — 10-20 in central Au-Au collisions) implies a much enhanced charmomium 
yield as compared to early (hard) production coupled with nuclear and QGP suppression mechanisms. However, 
at SPS energies, expected plasma lifetimes of ^ 1-2 fni/c, together with initial temperatures below ^ 250 MeV, 
may not suffice for a (close to) complete destruction of primordially produced of J/^" mesons. At the same time, 
a significant increase in open-charm production is not easily justified. This led us to propose a two-component 
model for charmonium production in heavy-ion collisions. It incorporates both a primordial yield (subject to 
subsequent dissociation) as well as a thermal contribution from statistical recombination of c and c quarks at the 
hadronization transition. Both components are evaluated within a common thermal fireball framework which is 
consistent with measured hadro-chemistry, ex pan sion dynamics, and has also been successfully employed to describe 
electromagnetic observables at SPS energies Another important feature is that we refrain from invoking any 

"anomalous" open-charm enhancement. 

In the present article, we expand upon our previous analysis in several respects: First, and most importantly, we 
give a detailed account of centrality dependencies for J/^ yields in both S-U and Pb-Pb systems at SPS, as well as 
predictions for Au-Au at full RHIC energy. The relevance of hadronic dissociation is assessed, especially in the context 
of the \E''/^' ratio measured at SPS. Higher charmonium states relevant for feeddown contributions to the J/4' are 
included on an equal footing. Effects of incomplete charm-quark thermalization, which affect thermal recombination, 
are incorporated on a phenomenological basis. We also lay out further steps to investigate problems that will not be 
satisfactorily addressed here. 

The paper is organized as follows: In Sect. II, we recall basic elements of nuclear absorption including normalization 



issues when comparing to NA38 and NA50 data. In Sect. IH, we present our calculations for the charmonium 
dissociation processes, i.e., parton-induced destruction in a QGP as well as inelastic hadronic J/^* interactions using 
effective Lagrangians. Pertinent survival probabilites are obtained from a convolution of the destruction rates over 
the space-time history of a collision modeled within an expanding fireball. Sect. ^ is devoted to the description of 
statistical charmonium formation at the hadronization transition including exact charm conservation constraints. In 
Sect. 0, we combine the two sources of charmonia to investigate the J/^ centrality dependence and 4*'/^ ratio at 



SPS for both S'(200 AGeV)-U and F6(158 AGeV)-P& systems. In Sect. 0, we give predictions for the ratio Nj/^/N^c 
as a function of centrality at RHI C, an d discuss the excitation function of this ratio from SPS to RHIC energies. We 
summarize and conclude in Sect. VIII| , and indicate directions for future improvements. 



II. NUCLEAR ABSORPTION OF CHARMONIA 

A suppressed abundance of J/'i> mesons relative to expectations from N-N collisions was first observed in proton- 
nucleus (p-A) and light-ion induced collisions. In these reactions, one does not expect noticeable effects from secondary 
particles. Indeed, rather extensive experimental studies |^, ||, |6| have established that the suppression can be 
understood as the absorption of a "pre-resonance" cc state in the (normal) nuclear medium of target and projectile 
nuclei. This in particular accounts for the observation that different charmonium states such as J/^* and 5"' follow 
very similar absorption patterns, despite their different bound state properties {e.g., binding energy or size). 

Let us briefly summarize the main elements of this "nuclear absorption" . In the Glauber model of A-B collisions, 
the probability for the pre-resonance state to be absorbed on its way through the nucleus is given by pH 

Snuc{b, (Tnuc) = „ ^ [ (f s dz dz' pa{s, z) pB{b - s, z') 



TAB{b) 

X exp|-(A-l)^ dzAPA(s, 2;a)(t««c| X exp |-(B - l)y" dzsPsib ~ s, ZB)(Jnuc^ (1) 

where pA (s , z) describes nuclear density profiles (taken from Ref. ) , s the position of the cc production point in 
the transverse plane, and b the impact parameter. The coordinates z and z' specify the positions within nucleus A 



3 



and B, respectively, along the collision axis. Tab (6) is the usual nuclear overlap function of the two colliding nuclei. 
The pre-resonance "absorption" cross section (Tnuc is treated as a free parameter being adjusted to experimental data. 

In the NA38/NA50 experiments, from which all data shown below are taken, the main measure for the centrality 
of a nuclear collisions is given by the transverse energy Et detected in the calorimeter. The impact parameter b is 
commonly related to this observable through the wounded nucleon model, which, upon inclusion of fluctuations in 
Et at fixed 5, provides a correlation function 

•p ^^ 1 / [ET-qNpart{b)l_\ 

Here, g is a proportionality factor between Et and the number of participant, ET(b) = qNpart, which depends on the 
specific experimental settings, and a is an empirical parameter characterizing the magnitude of the fiuctuations. The 
absorption of charmonium in nuclear matter as a function of Et then follows as 

Bf^f^dj/xj, _ Bf^f^Cr"/^ J (PbVAB{ET,b)Snuc{b,(Jnuc)TAB{b) 

<JDY ~ CjPPy jdPbVAB{ETMTAB{b) ^' 

where B^^ is the branching ratio for J/^* — > fi^ fi^ . The pre-factor cr^'j^^^/a^y is not accurately known, and depends 
on the energy of the collision. 

At SPS energies, we will concentrate on the two collision systems S'(200 AGeV)-C/ and P6(158AGeV)-P6. For q 
and a we use the values reported by the experiments, i.e., q = 0.275 (0.72) and a = 1.27 (1.56) for Pb-Pb {S-U). 
Together with an absorption cross section of anuc = 6.4 mb and normalization factors of 52.8 for Pb-Pb and 48.2 
for S-U, our results coincide with nuclear absorption studies of NA38 and NA50, cf. Fig. |l|. Whereas the S-U data 
are well accounted for, the Pb-Pb system departs from nuclear absorption starting at Et — 40 GeV, constituting 



"anomalous" J/'^ suppression l23, p4, 25, Era 




FIG. 1: Nuclear absorption calculations for the S'(200 AGeV)-U (left panel) and P6(158 AGeV)-P6 (right panel) systems at 
the CERN-SPS, compared to data from NA38 |^ and NA50 |2^], respectively. The normalization factor B ^^j^a'^l^ / has 
been fixed at 48.2 (52.8) for the S-U (Pb-Pb) system. 



III. J/* DISSOCIATION IN HEAVY-ION COLLISIONS 



The deviation of the J/'i> yield from nuclear absorption systematics in Pb{158 AGeV)-P6 has triggered extensive 
theoretical analyses. The responsible underlying mechanisms are still a matter of debate, ranging from destruction 
in a QGP and/or hadron gas to thermal production at Tc. As stressed above, this work attempts a comprehensive 
treatment of all three of these aspects within a minimal set of assumptions, which, in particular, is based on a thermal 
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description of the collision dynamics including Quark-Gluon Plasma formation if the initial conditions are energetic 
enough. In this section, wc compute J /'^ survival probabilities for both the plasma and hadronic phases of heavy-ion 
reactions. 



A. Quark-Gluon Plasma 



J /'^ dissociation in a thermalized QGP has been studied in both static and dynamical frameworks. Within the 
former, one typically evaluates the screening of the heavy quark potential by color charges, whereas the latter involves 
inelastic parton collisions using, e.g., the QCD analogue of photo-dissociation, g + J /'^ — > c -I- c [^ . 

Here we adopt a dynamical approach, accounting however for reduced charmonium binding energies as extracted 
from a Schrodinger equation for cc bound states in a screened heavy-quark potential [ p8| . Under conditions relevant 
to the initial stage of heavy-ion collisions, the binding energy of the J/5' is strongly reduced with respect to its 
vacuum value, which is also borne out of recent lattice gauge calculations ]29[| . Color-screening is characterized by the 
electric Debye mass /iu which wc estimate to leading order in perturbation theory, /iu ~ gT. Since, strictly speaking, 
perturbation theory is not really applicable under the moderate plasma temperatures expected even at RHIC energies, 
we regard the strong coupling constant as an effective parameter to be adjusted to the J/^* data at SPS. E.g., with 
a typical g ~ 1.7, the J/-^ binding energy at T = 170 MeV is Edtss = 250 MeV decreasing to Ediss = 100 MeV at 
T — 230 MeV and vanishing around T ~ 360 MeV. For such small binding energies, the photo-dissociation process 
becomes inefficient due to unfavorable break-up kinematics. For a loosely bound charmonium state, inelastic parton 
scattering, g{q, q) -f J/"^ g{q, q) + c + c, turns out to be a more important mechanism |l9[| . The respective cross 
sections are evaluated in quasifree approximation using leading-order QCD for gc — )■ gc {qc — > qc) ]30[ |. The thermal 
dissociation rate is then obtained via 



E 



rfik 

—^nk,T)adUs) (4) 



(27r)3 



denotes the minimal on-shell momentum of a quark or gluon from the heat bath necessary to dissolve an 
in-medium charmonium bound state into a free cc pair. We include thermal quasiparticle masses for light quarks and 
gluons 

TO„ = mn H m„ = — — . (5) 



6 ' " " 6 
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This formalism straightforwardly applies to ^' and x states as well. From the Schrodinger equation, it follows |28| that 
their binding energy vanishes at or even below Tc, which is also found in recent lattice calculations Therefore, the 
pertinent QGP dissociation rate is given by Eq. (jj) with Ediss = 0. The picture we have in mind here is similar in spirit 
to that underlying nuclear absorption: a pre-resonance state, characterized by a correlation induced in the primordial 
A'^-A'^ collision {e.g., comoving c and c quarks), is only destroyed if it actually undergoes an inelastic interaction with 
the surrounding matter. 

The resulting charmonium dissociation times, Tdiss — ^diss^ shown in Fig. |^ and, for the J/^*, are compared to 
photo-dissociation without medium effects in the bound state energy (as has been employed in the literature before). 
For temperatures relevant at SPS, the quasifree process (full line) is more efficient in dissolving J/5''s than photo- 
dissociation (dotted line). At higher temperatures this tendency is reversed, mostly due to an increasing Debye mass 
which suppresses the <-channel gluon-exchange graphs for g{q, q) + J/^ — > g{q, g) -I- c -I- c. Close to Tc, the vanishing 
binding energies for and x entail lifetimes which are about a factor of three below the one for J/\t'-mesons. 



B. Hadron Gas 



Subsequent to quark-gluon plasma dissociation, when the system converts to the hadronic phase, charmonia undergo 
further suppression due to interactions with surrounding hadrons. There are a number of theoretical models for 
J/^* -I- hadron processes |3^, |^ |3^, ^ whose results span an appreciable magnitude. Calculations 

involving excited charmonia, such as ^E"' and Xc (which contribute to the J/^* yield via electromagnetic feeddown), are 
scarce. Consequently, the impact of inelastic hadronic scattering is not very well under control. Most calculations, 
including recent ones based on quark exchange models ||3^ or on S'[/(4)-symmetric effective lagrangians |39|, seem 
to indicate a rather moderate impact of the hadronic medium on the J/^!. 
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T [GeV] 

FIG. 2: Dissociation times for a J/'^ , x) in a QGP as a function of temperature. The full (dashed) curve corresponds to 
the leading-order QCD process for quasifree g,q + c ^ g,q + c scattering with in-medium J/'i/ , x) bound state energies. 
The dotted curve results from photo-dissociation, gJ/^ — > cc, assuming the vacuum dissociation energy. 

1. Absorption on light hadrons 

As a baseline calculation, we here reproduce results obtained in Refs. ^ within a S'C/(4) effective theory. The 
starting point is a S'[/(4)-flavor symmetric effective lagrangian formulated with 4-by-4 pseudoscalar and vector meson 
matrices. Although the SU (4) symmetry is strongly broken by the c-quark mass, the hope is that symmetry-breaking 
effects are largely accounted for by the hadronic mass matrix. Along these lines we compute (inelastic) interactions 
of the J/* with pions (7r + J/* D + D*,D + D*) and rhos {p+J/^ ^ D + D, p + J/'i' D* + D*) which are the 
most abundant mesons in the medium. We employ coupling constants as calibrated by Haglin and Gale in Ref. ]39| ] 
based on the p tttt decay to fix the gauge coupling. The effective hadronic theory is supplemented by vertex form 
factors to simulate finite-size effects. This violates gauge invariance, which, however, can be restored by introducing 
appropriate counter terms |39|. We have checked for various processes {e.g., tt -|- J/^ D + D*) that these extra 
terms induce only quantitatively minor modifications, which we therefore neglected. The corresponding J/^* lifetime 
in a TT-p gas is displayed in Fig. using (covariant) monopole form factors with cutoff A = 1 GeV. 




T [GeV] 

FIG. 3: Hadronic lifetimes for J/^ (solid line), ^' (dashed line) and x (dot-dashed line) based on a S'C/(4) effective lagrangian, 
including interactions with pions and rhos using a form factor value of A = 1 GeV. The curves for the (x) axe deduced from 
the J/* one via cross section scaling by (r'*'/''J/*)^ ((^x/''J/*)^ for the x)- 

As mentioned above, and x dissociation rates are even more difficult to assess within hadronic model frameworks. 
To obtain an estimate that can be used in our calculations below, we assume a geometric scaling of the calculated J/^if 
cross sections with the squared ratio of the respective charmonium radii (?'*'/'".//* )^ ^^I'i (^x/^.//*)^' dashed 
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and dot-dashed lines in Fig. ||. It turns out that these estimates agree reasonably well with explicit calculations of 
the corresponding processes within the constituent quark model approach of Ref . |]38| . 



2. Anomalous Processes 



We furthermore investigated the role of so-called anomalous processes such as 7r + J/^ — > ric + p, 7r + J/^ f]c + bi, 
and p + J/"^ rjc + TT, suggested in Ref. |39|. In there, the relevant hadronic couphng constants, e.g., gj/^uria^ have 
been estimated applying the vector dominance model (VDM) to the radiative decay J/^ ^ 'yrjc- We believe, however, 
that this procedure leads to a significant overprediction of the hadronic coupling, for the following reason. For vertices 
carrying identical quantum number structure, namely J I'^Lorj and J /"^wq' , both hadronic and corresponding radiative 



decay information is available from experiment. One finds |40 



■ X 10 
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r( J/* ^ uji^)/Ttot = 1.6 X 10-3, r( J/* ^ uJ7^')/Ttot - 1.7 , , _ 
r( J/* ^ lv)lTtot = 8.6 X 10-^ r( J/* iri')lTtot = 4.3 x IQ-^ , 

i.e., the hadronic branching ratios are comparable or even below the radiative ones. This is in marked contradiction 
to VDM within which the latter are suppressed by a factor iej g^)"^ which is much smaller than the moderate increase 
in phase space due to the final state decay momenta (also note that VDM is more strongly violated with increasing 
mass of the pseudoscalar meson). Thus it appears that VDM cannot be applied, and that an accordingly reduced 
Sj/'i'ujije coupling renders the pertinent ^-channel w exchange processes in tt -|- J/\l/ ^ rjc + p and tt + J/'^ ^ ijc + bi 
negligible. We therefore decided not to include anomalous processes in our analysis. 

C. J/'i survival probability in a thermal expansion scenario 

The final number of primordial J/^I^'s remaining after plasma and hadronic phases requires the convolution of the 
dissociation rates derived in the previous two sections over the space-time history of a given heavy-ion collision. To 
this end, we model this evolution by a schematic thermal fireball expansion ]4l| ] , which, however, incorporates essential 
features of hydrodynamical calculations. Let us briefly recall the main elements. 

Equilibration of the system is assumed at a formation time tq, after which isentropic expansion proceeds at flxed 
entropy per baryon, S/Nb = s/ub, where the total entropy S = sVfb and net baryon number Nb — tlbVfb are 
related to the pertinent densities via the time-dependent 3- volume Vfb ■ The latter is modeled in cylindrical symmetry 
as 

Vfb{t) = 2(zo + v^T + ^a,T^)Tr{ro + (6) 

(the overall prefactor of 2 accounts for two fireballs which enable to cover about 4 units of rapidity), rg denotes the 
initial transverse overlap of the two colliding nuclei at given impact parameter &, whereas the expansion parameters 
{vz,az,a±} are adjusted in line with hydro-calculations to reproduce observed flow velocities in connection with 
(thermal) freezeout times of Tfo — 10—14 fm/c. The parameter zq is equivalent to the formation time tq (in the 
Bjorken limit zq ~ tqAj/), specifying the initial conditions of the evolution. The temperature of the expanding matter 
is inferred from the entropy density s(r) = S'/Vfb(''") in either hadronic or QGP phase (calculated from an ideal 
resonance gas for the former and using massive quasiparticles according to Eqs. (||) for the latter). If s(t) lies in 
between the entropy densities and s^'^ for hadron gas and QGP, a standard mixed phase construction is 
employed at the critical temperature T^, 

/^f + (7) 



Vfb it) 



which determines the volume partitions / and 1 — / for hadronic and quark-gluon matter, respectively. 

The collision-energy dependence of the underlying parameters from SPS to RHIC has been flxed as follows. The 
total entropy at given collision energy and centrality is obtained from the speciflc entropy S/Nb ranging from 26 to 
250 from SPS to RHIC according to standard chemical freezeout analyses [|[ ^ (with effectively conserved numbers 
of pions, kaons etc., thereafter). The critical temperature is assumed to increase smoothly from 170 MeV at SPS 
to 180 MeV at RHIC (with ^b decreasing from 260 to 27 MeV). In addition, we account for 10-20 % variations in 
S/Nb with centrality following particle production systematics reported by NA49 j4^ and PHOBOS |Q. Finally, 
the formation time tq (i.e., zo) is continuously decreased from 1 to 1/3 fm/c for ^/snn from 17 to 200 GeV. With 
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the volume expansion as given by Eq. (^), this resuhs in QGP phase durations in central collisions between 1 and 
~ 3 fni/c followed by mixed phases lasting 3-4 fm/c. 

We are now in position to calculate the desired survival probabilities. The evolution equation for the number of 
each charmonium species i (i = J/vl/, VE*', x) iii the system at time r reads 



dr 



— — 



The destruction rate T^^^^ is specified according to the phase and temperature T(t) of the system, 



QG ' 

d^ss = { /rk + (i-/)rbG 



T>T, 
T<T, 



(8) 



(9) 



(with / given by eq. (|^), and and from Sects. Ill A and IIIB). Eq. is readily integrated to obtain the 
survival probability at time r during the collision, 



(10) 



The final survival probability, Sqg+h, relevant for experimental observables is simply given by the value of Sgc+Hi''') 
at the moment of thermal freezeout, t/o, where all hadrons cease to interact. Note that '5qq_|_^ depends on the impact 
parameter b of the collision through the space-time evolution of the system. 

Fig. ^ shows J/^" and survival probabilities at full SPS and RHIC energies (central collisions, Npart = 360) as a 
function of time with a value of 5 = 1.7 for the strong coupling constant (which provides reasonable agreement with 
SPS data to be discussed below). Most of the suppression originates from the plasma and mixed phase, with hadronic 




Au-Au Vs^j^=200 GeV 



Au-Au Vsj^^=200 GeV 



J/l*, Pb(158 AGeV)-Pb 
'i", Pb(158 AGeV)-Pb 



5 10 

time [fm/c] 
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FIG. 4: Survival probability of J/* (dot-dashed line: SPS, solid line: RHIC) and *' (dotted line: SPS, dashed line: RHIC) 
for central collisio ns (Np art ~ 360 ) as a function of time. The curves are obtained upon integration of the dissociation rates 
calculated in Sect. [HA and [II B| . 



effects playing little role. The destruction of ^P' states (dotted line) under SPS conditions is far from complete (with 
a fraction of about 1/3 remaining at the end of the collision); this may in fact point at a lack of our understan ding 
of the hadronic 5*' interactions, since (i) it will cause problems in describing the measured 5''/4' ratio (see Sect. |VB| 
below), and (ii) it appears to be incompatible with recent lattice calculations which exhibit a dissolution of the 
in the hadronic phase (although no statement about timescales could be made). As expected, the overall dissociation 
is much stronger at RHIC (full line: J/^t dashed line: 5*'), although 25% of primordial J/^'s endure the evolution. 

For a combination of the (suppressed) primordial production with the statistical one (to be discussed in the next 
section) it is necessary to convert the survival probabilities into absolute yields. At a given impact parameter b, the 
number of J/^' mesons initially produced in hard parton-parton collisions is 



N°/^{b)^ail^ABTAB{b) . 



(11) 
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There are no data on charmonium production in p-p in the RHIC energy regime (to date, the highest cms energies 
available lie around ^/s'nn — 40 GeV). Thus we have to rely on extrapolations. The so-called Schuler parameteriza- 
tion Q - a phenomenological fit based on low-energy systematics - leads to fTp/*(-\/s = 200 GeV) = 1.0 /^b. On the 
other hand, the simple ansatz CTp/* = fcTpp with / ~ 0.025 Q together with PYTHIA extrapolations of open-charm 

production (see also next section), which also reproduces low-energy data, gives Opl"^ {y/s = 200 GeV) = 8.75 ^b. 
This obviously introduces appreciable uncertainties in establishing the J/^* yield arising from primordial production 
mechanisms. 

In addition, higher charmonium states (in particular 5"' and various x states) contribute significantly to the measured 
J/vP abundance via strong and electromagnetic decays (feeddown). Since these resonances have different hadronic 
and plasma cross sections, the distinction between prompt and secondary J/^'s has to be made. Consequently, the 
total number of J/^'s stemming from primordial charmonium states after nuclear, plasma and hadronic suppression, 
becomes 

Nff^ih) = aJJ^ABTABib)S„uc [0.65^/,% + O.OSS^g+h + O-^^S^g+h] ■ (12) 



IV. STATISTICAL PRODUCTION OF J/* 

Besides direct J/^E* production, another source of J/^I^'s has been attributed to thermal (or statistical) production in 
several variants. We here adopt the approach put forward in Refs. [|2], , where statistical coalescence of primordial c 
and c quarks is assumed to predominantly occur at the hadronization transition. The relative abundances of charmed 
particles are then determined by hadronic thermal weights at the critical temperature Tc, whereas the absolute number 
is specified by the amount of charm quarks created in the early (hard) TV- TV collisions. 

Starting from the ideal gas expression for the thermal number density of particle species i, 



2^ 



p dp 



exp 



m, - pLi 



T 



± 1 



(13) 



{di. spin-isospin degeneracy), the total number densities of open and hidden charm particles follow as 



(14) 



rihid 



(15) 



The pertinent chemical potentials encode conserved charges [net baryon, strangeness and charm numbers), and 
depend on the collision energy of the system (in accord with the hadro-chemistry) . At we match the densities, 
Eqs. ( p^ ) and (|l^), to the available number of open charm pairs, N^c, by introducing an effective fugacity factor 
7c ~ 7c for both charm and anticharm quarks. Since A'cg is typically a small number (e.^., ^0.2 for central Pb-Pb at 
SPS), exact charm conservation within the canonical ensemble formalism Q is mandatory, leading to 



Ncc = -IcnopVH—, ^^7^ + ^c'^MciVh ■ (16) 

2 loi^cnopVH) 



Here, Vh denotes the fireball volume from the previous section, cf. Eq. (^, at the moment when hadronization is 
complete. The actual value of Ncs depends on both collision energy and impact parameter. The former dependence 
is inferred from PYTHIA computations upscaled by an empirical if-factor of around 5 extracted from a best fit 
to existing pN and irN data |5^. As discussed in the previous section, the extrapolation into the RHIC regime bears 
appreciable uncertainty. E.g., next-to-leading order pQCD calculations for full RHIC energy give ~ 350 ^b, 
to be compared to app ~ 570 ub from the PYTHIA estimate. Both extrapolations are in line with recent indirect 
measurements from PHENIX in Au-Au at y^SNN = 130 GeV, where single-electron pt-spectra have been used to 
infer ct^^ — 380 ± 60(stat) ± 200(syst) /xb, to be compared to app ~ 320 /ib within the PYTHIA extrapolation. 

By construction, statistical charmonium production is only active for c and c that emerge from a deconfined 
environment prior to recombination. Therefore, in peripheral collisions where the initial volume is only partially filled 
with the QGP phase, only a fraction of the cc pairs is available for coalescence. Accordingly, the hadronization volume 
which enters Eq. ( [l6|) arises from the initial volume in the plasma phase after expansion. 
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The thermal equiUbrium (but chemical ofF-equilibrium) J/^t yield takes the form (J/^) — 7^ Vh nj/^i,. The former 
expression is valid if (anti-) charm quarks are kinetically equilibrated, i.e., the momentum distribution of c and c 
quarks is thermal which is questionable under SPS and even RHIC conditions. We therefore implement the following 
correction: we introduce a thermalization time Teq for c and c quarks, approximated as r^q = 1/na, where n is the 
total density of quark- and gluon-quasiparticles in the system, and a is the elastic scattering cross section for the 
processes g{q, q) + c(c) —^ g{q, q) + c(c). Within a relaxation time approach, the relative reduction TZ in thermal J/^ 
formation is then estimated as 



n 







1 — exp 















(17) 



where th is the time at which hadronization is completed. Our expression for the number of J/^''s produced at the 
hadronization transition by coalescence of a c and c quark is thus modified according to 

(J/^) = 7f Vh nj/^n . (18) 

We have checked that our results do not change significantly with the definition adopted for th. The latter is not a 
sharply defined quantity since the mixed phase lasts for a few fm/c. However, smaller values of th {e.g., if taken in 
the middle of the mixed phase) lead to smaller volumes and larger implying an increase in thermal J/^* production 
which, in turn, is (partially) compensated by a less degree of thermalization through a smaller value of TZ. 

The charmonia statistically produced at the hadronization transition are still subject to reinteractions in the 
hadronic phase, so that their final contribution to the observed J/4' yield is accounted for via 



X 



J/* + BR{X ^ J/*) nx 5f 



n , (19) 



where the summation is carried over the charmonium states X with finite decay branching into J/4''s (feeddown). 

V. TWO-COMPONENT MODEL AT THE SPS 

By combining the two sources of charmonium as discussed in the previous sections, we now turn to applications 
to heavy-ion collisions at SPS energy (v^jvat ~ ^'''-^ GeV). We first address the observable that has drawn the most 
attention, i.e., the centrality dependence of the J/^* yield, but also investigate the 4''/\E' ratio. 

A. Centrality dependence of ,J/^ production 

One of the important findings of the NA38 and NA50 experiments at CERN is that the J/^* yield in p-p and 
p-A collisions with light and heavy targets is well explained by hard production coupled with nuclear absorption, 
cf. Sect. II. Thus, any attempt to describe the J/^ yield in heavy-ion collisions must reproduce this feature. In 
particular, for very peripheral collisions involving only a few participant nucleons, the J /'^ suppression pattern should 
coincide with nuclear absorption. 

In our two-component model, the total number of observed J/4''s is the sum of direct and statistical production 
according to Eqs. ( p^ and (|l^), respectively, 

Nj/^ = Nff^+Nff^ . (20) 

Both terms on the right-hand side depend on centrality (via the impact parameter h) and collision-energy. For very 
peripheral collisions in the SPS regime, the initial conditions are not energetic enough to induce a transition into the 
QGP. Therefore, the J/'^ source of statistical recombination at Tc is absent; at the same time the hadronic interactions 
are not frequent enough to induce sizable effects. Thus, for large impact parameters our approach is consistent with 
the NA38/NA50 results. 

For a detailed comparison with data, we need to evaluate the ratio B^^^a^ 1^ j commonly displayed by 
NA38/NA50 as a function of transverse energy deposited in their calorimeter. Following the treatment outlined 
in Sect. II, we convolute our theoretical results (which are obtained as a function of impact parameter h) with the 
probability distribution V{ET,b). This amounts to replacing the numerator in Eq. (|^) by the expression 

(Et) = B^^aPP^ I d% Vab{Et, b) [s„ucSqg+h + Nf^^ / {aPP^ABTAB{b))\ TAeib) . (21) 
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Our results are compared to NA38/NA50 data in Fig. | for both the 5(200 AGeV)-U (left panel) and P6(158 AGeV)- 
Pb (right panel) system. At all centralities, direct production (dashed line) prevails over the thermal component 




FIG. 5: Results on J/^' nroduction within the two-component model as a function of centrality compared to NA38 g (left 
panel) and NA50 data ^ for 5(200 AGeV)-U and P6(158 AGeV)-Pb, respectively. Dotted lines: direct production with 
nuclear absorption alone; dashed lines: direct production subject to nuclear, QGP and hadronic absorption; dot-dashed lines: 
statistical (thermal) production from a hadronizing QGP including hadronic dissociation; solid lines: combined direct and 
statistical yield (sum of dashed and dot-dashed lines). 

(dot-dashed Hne). The latter sets in once a QGP starts forming, which, in turn, requires a stronger QGP suppression 
of the direct component than without the thermal contribution. The adjustment of the only free parameter (strong 
coupling constant g = 1.7) to the most central Pb-Pb data allows for a satisfactory reproduction of the centrality 
dependence for this system. For S-U collisions, the results are somewhat on the low side. Note that in our approach 
the "drop" in the Pb-Pb data around Et — 40 GeV is a combination of a rather strong QGP suppression coupled 
with the onset of thermal production. 

In its present form, our model does not capture the appearance of the "second drop" in the data for the most 
central Pb-Pb collisions at Et > 100 GeV. In fact, the maximum transverse energy in our description is at E™'^^ = 
E-rib = 0) = 100 GeV, well below the experimental limit which extends up to Et — 125 GeV. It has been suggested 
that these features are associated with transverse energy fluctuations H^, |5^ and/or trigger energy losses and 
are thus not necessarily related to a shortcoming in an underlying (microscopic) model description of J/^! production 
(suppression) . 

Let us first address the Et fluctuations. From the minimum bias (MB) event distribution of transverse energy, 
dN/dET, as measured in the NA50 apparatus |Q one finds a rapid falloff beyond Et = 100 GeV, the so-called 
"knee" of the distribution. The tail of the latter, which reaches beyond Et ~ 100 GeV, is associated with fluctuations 
in transverse energy at fixed geometry for 6 = 0. Events which fluctuate beyond Et > 100 GeV contain a larger 
initial entropy and consequently a hotter initial temperature and a longer plasma phase which implies additional J/^> 
suppression (charm and J/^> production being a hard process are not coupled to fluctuations in the "soft" sector). To 
account for this phenomenon we replace in our calculations the total entropy at fixed impact parameter, Stotib), 
by 

Stotib) ^ Stotib)^-^ . (22) 

This does not affect our results for Et < ET{b = 0) ~ 100 GeV, but beyond the knee of the distribution, Stot is 
enhanced by the factor Et / ET{b). The effect of this modification is shown in the left panel of Fig. |[ where we 
compare our calculations to data on J/^ production normalized to the minimum bias £^T-distribution^ . Obviously, 



^ This way of normalizing the data has the advantage of the much better statistics for dN/dEx as compared to the Drell-Yan sample. On 
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the description of the observed turnover for Et > 100 GeV is improved by inclusion of £'T-fluctuations (cf. dashed 
versus dot-dashed curve), but does not suffice to quantitatively explain the data. Note that this effect relies on 
additional J/^! suppression in both the direct component (due to stronger suppression) and the thermal component 
(due to a larger hadronization volume weakening the enhancement induced by the canonical-ensemble treatment). 




FIG. 6: Left panel: J/^ over Minimum Bias [MB) ratio compared to calculations within the two-component model (dot-dashed 
curve) with additional inclusion of effects due to Et fluctuations (dashed curve) and trigger energy loss (solid curve). Right 
panel: impact of a model for trigger energy loss on the Drell-Yan (DY) over MB ratio. All data points are from NA50 for 
P6(158 AGeV)-P6 ||^. 

However, as has been suggested by Capella et al. ||5^, there might be an additional feature in the large i5r-region 
unrelated to J/^' physics, which can be gleaned from the Drell-Yan (DY) over minimum bias data, cf. right panel 
of Fig. H The theoretical ratio DY/MB (dotted line) computed by NA50 flattens out at large values of Et > 100 
GeV, whereas the data seem to indicate a slight turnover |p^ . The argument to explain this turnover (which 
equally applies to the J/^ event sample) is that for DY-events the hadronic transverse energy deposited in the 
calorimeter is slightly reduced as compared to the corresponding Afi?-events due to triggering on the DY (or J/^*) 
pair (2.9 GeV< Mdy <4.5 GeV). As elaborated in Ref. ||5^, a rough estimate of this effect is obtained by rescaling 
the amount of transverse energy in the J/"^ and DY events according to 

i?H6)->£;^(6)%i^. (23) 

When incorporated into the Ex-b correlation, Eq. (||) , this rather small loss in Et entails a surprisingly large drop of 
about ~20-30% in the tail of the (DY/MB) distribution, cf. Fig. ^ (solid curve in the right panel). 

Returning to the {J/'^)/MB ratio (left panel of Fig. H), we see that the combined effect of Et fluctuations and 
trigger energy loss, implemented within our microscopic model, gives a satisfactory description of the experimental 
findings. This also holds true for the more common representation of the data via the J jDY ratio, cf. Fig. |^. Note 
that in this case the trigger energy-loss correction only applies to the "MB" data sets, which have been extracted 
using a theoretical expression for the MB/ DY ratio according to 

DY ) MB analy s^s \ M B J \ DY J ^ ' 



B. *7* ratio 



In p-p collisions in the SPS energy regime, the ratio of produced ^P' to J/^ mesons amounts to a value of about 
15%, which persists for p-A collisions as nuclear absorption affects both charmonium states in practically the same 



the other hand, since Et production is governed by soft physics (which essentially scales with the number of participants rather than 
the number of A'^-A' collisions), the characteristic features of J/^ suppression are not readily discernible. 
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FIG. 7: Results of the two-component model without (dot-dashed line) and with additional inclusion of transverse energy 
fluctuations (dashed line) and trigger energy loss (full line), for the centrality dependence of the _B^,,cr''/*/(T^^ ratio in P6(158 
AGeV)-Pfe colhsions. 



way, cf. Sect. II. A marked deviation from this behavior has been observed in S'(200 AGeV)-C/ and P6(158 AGeV)-P6 
collisions, with an onset at rather low centrahties. Remarkably, for central collisions, the ^f'/^ ratio does not go to 
zero, but rather levels off at a value of around 3-4%. This is contrary to the naive expectation that ^P' mesons, due 
to their much smaller binding energy than the J/^t states, are significantly more suppressed. In Ref. ||5^ it has been 
suggested that, under the premise that most states are dissolved in the QGP, their abundance is regenerated in the 
hadronic phase from remaining ^> states as a consequence of chiral symmetry restoration, via the process ^'-I-tt ^ 4''. 
The interaction was assumed to be mediated through cr(500) meson exchange, the mass of the latter approaching the 
pion mass thereby substantially enhancing 4"' formation. From another point of view, the fact that the value of 4% 
reflects the thermal ratio of \E''/\1/ at a temperature of T = 170 MeV (with little latitude), has been put forward in 
ref. as evidence for statistical charmonium production at the hadronization transition. 

Within our two-component model as laid out above, the ^f'/^" ratio follows without further assumptions. The 
results for both S-U and Pb-Pb systems are compared to NA38/50 data in Fig. ||. The discrepancy with experiment 



□ Pb-Pb 96 




FIG. 8: Our calculations for the centrality dependence of compared to data from NA38/50. The solid and dashed line are 

for 5(200 AGeV)-U and Pfe(158 AGeV)-P6, respectively, using hadronic 'I'' dissociation cross sections obtained by geometric 
scaling of the J/'i! one, cf. Fig. ^ The dashed and dotted lines are the corresponding results when artificially increasing the 
'5' cross sections by another factor of 5. 
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is rather significant, especially for 5(200 AGeV)-J7. Since in the latter reaction QGP effects are not expected to play 
a pronounced role, it seems that the deficiency in our description has to be assigned to the hadronic phase, i.e., an 
underestimation of the hadronic cross sections for the Indeed, an artificial increase of this quantity by, say, a 
factor of 5 clearly improves the agreement with the data. We have checked that such an increase in the hadronic 
cross sections has negligible impact on the "^/DY ratio as plotted in Fig. |^ (the contributes maximally 8% to the 
observed J/vP yield). 

As mentioned before, recent lattice calculations indicate a dissolution of the ^f' in a static environment well below 
the phase transition temperature, due to the lowering of the DD continuum threshold below the in-medium 5"' mass. 
In a hadronic model framework, this can be implemented by an in-medium reduction of the D-meson masses, which 
has been motivated in Ref. p9[ by chiral restoration arguments inducing a lowering of the light-quark mass within 
the cq and cq states. We have investigated this possibility within the hadronic approach employed here, and found a 
strong sensitivity to the detailed modeling of the light-quark related portion of the Z3-meson masses. In fact, if the 
lifetimes of the charmonium states become comparable to duration of the fireball expansion, one needs to account for 
the reverse reaction of charmonium formation (as required by detailed balance), which is beyond the scope of this 
paper^ . 

A more controlled way to assess hadronic medium effects in charmonium dissociation should be provided by con- 
stituent quark (exchange-) models incorporating both phenomenological confinement potentials as well as properties 
of chiral symmetry breaking ^ . This will be addressed elsewhere . 



VI. EXCITATION FUNCTION AND PREDICTIONS FOR RHIC 



An essential part of the experimental program at RHIC is again on (penetrating) electromagnetic probes. The 
PHENIX detector will provide accurate dilepton data via both the (forward) muon arms as well as electron identifica- 
tion in the central region. The results on charmonium should allow for stringent constraints on models. At full RHIC 
energy, standard extrapolations predict an open-charm production that is about two orders of magnitude larger than 
in the SPS regime, entailing a substantial increase in the statistical recombination mechanism for charmonia. At the 
same time, direct (hard) charmonium production, albeit also enhanced by presumably a similar factor as open charm, 
ought to be more strongly suppressed due to longer and initially hotter QGP phases. 

A quantitative comparison between SPS and RHIC within our two-component model is performed in Fig. |p where 
the final (observed) number of J/'^'s, normalized to the number J/\l/'s remaining after nuclear absorption, N^J^, is 
displayed for central collisions as a function of the fireball evolution time. 
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FIG. 9: Time dependence of the ratio Nj/ig/NjJ"^ at SPS (dashed line) and RHIC (full line) for central collisions with 
Npart = 360, where NjJ"^ is the number of J/5''s remaining after nuclear absorption. The respective fractions of direct (/dir) 
and thermal (fth) yields are indicated by the arrows. 



^ Note that with an increase by a factor of 5 for dissociation rate over the results shown in Fig. ^ (as applied in Fig. g) the ^' lifetimes 
in the vicinity of Tc are indeed close to the expansion time of the hadronic phase. 
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The freezeout value for this ratio increases by about 50% going from SPS to RHIC (from 0.65 to about 1). More 
dramatically, the composition in terms of underlying sources is very different: whereas at SPS the (suppressed) direct 
yield dominates, J/^'-mesons at RHIC originate to ^80% from thermal production (being proportional to (TVcg)^). 
The upward jump of the two curves in Fig. ^ is located at the respective end of the mixed phase, th, where in our 
approximation all thermal production is concentrated. As elucidated in Sect. IV, the final results are not sensitive to 
the exact production time within the mixed phase. 

Since the phenomenological extrapolations for absolute numbers of primordial J/^! and charm-quark production 
up to RHIC energies, which are input parameters to our calculations, are beset with appreciable uncertainties, it is 
desirable to define a quantity which reduces this sensitivity. Therefore, we show in Fig. ^ the predictions of our 
two-component approach for centrality dependencies of the ratio Nj/^/Ncc, which, in anticipation of open-charm 
measurements at RHIC, also has the virtue of being composed of experimental observables'^. Whereas at SPS energies 
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FIG. 10: Comparison of the centrality dependence of the Nj/q,/Ncc ratio at SPS and RHIC full energy (full line). The direct 
(thermal) contributions are shown separately, respectively in the dashed (dot-dashed) lines. 

this quantity exhibits a monotonous decrease with increasing number of participants (left panel of Fig. ^), it saturates 
already for rather peripheral collisions at full RHIC energy. Note that the decrease of the thermal component for 
Npart > 100 at SPS is caused by canonical ensemble effects, while at RHIC statistical production, for the most 
part, proceeds in the grand-canonical limit entailing a smooth increase with centrality. Our approach thus clearly 
discriminates between standard J/\l/ suppression as opposed to thermal regeneration at full RHIC energy. 

It is therefore important (and experimentally feasible at RHIC) to map out the transition between the regimes of 
predominantly direct to thermal charmonium production, as has been first pointed out in Ref. p9| . In Fig. |ri| we 
present an updated prediction of the excitation function for Nj/^^, jNcc ratio"*. The ratio exhibits a nontrivial minimum 
structure around ^/snn — 40 GeV, which is a marked feature of the interplay between hard and thermal production 
(assuming no anomalies in open-charm production, cf. Ref. [ |64| ) . 

For practical purposes, it is also of interest to convert our results into an absolute number of J/5''s predicted for 
RHIC. Keeping in mind the uncertainties mentioned above, we find that the combined statistical and direct J/^" yield 
amounts to rapidity density dN/dy ~ 1x10^^ for central collisions at full RHIC energy. 



VII. COMPARISON TO OTHER WORKS 



The evaluation of J/'? dissolution in a QGP dates back to the late 70's / early 80's and has been discussed in many 
facets. In our approach, we have combined an in-medium reduced binding energy with parton-induced destruction. 



^ However, this ratio may still be sensitive to, e.g., shadowing effects: since the J/^ yield, being mostly thermal, goes as N^^ a modification 
of cc production due to nuclear shadowing does not cancel out in this ratio. 

As compared to our earlier results, we here also incorporate the corrections from feeddown of excited charmonia and from hadronic 
suppression. Within the uncertainty band quoted in Ref. |hg|, the results agree. 
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FIG. 11: Excitation function of the Nj/^, /Ncc. ratio (full line). The interplay between the direct contribution (dashed line) and 
the thermal component (dot-dashed line) results in a minimum in the excitation function around ^/snn — 40 GeV. 



which we have treated in quasifree approximation due to small (if any) binding energies of the various charmonium 
states. The latter feature is also the main origin of the rather strong suppression we find using rather moderate values 
of the strong coupling constant, as ~ 0.3, which, of course, is to be regarded as an effective parameter which we 
adjusted to describe the NA50 data. Due to a ~20% contribution from statistical J/\l/ production at SPS energies, 
our QGP suppression required to reproduce the NA50 data is stronger than in calculations based on the suppression 
effect alone. On the other hand, within the quasifree approximation, the discrimination between different charmonia is 
less pronounced than within "threshold" -type approaches where higher charmonium states are completely suppressed 
once the QGP temperature exceeds the relevant value. 

Concerning the evaluation of thermal (statistical) recombination, our approach parallels the one initiated in Refs. |lj, 
which is based on an open-charm abundance generated by hard production, available for recombination at Tc- 
The main difference to these analyses is that we do not invoke an enhancement of open-charm production beyond 
standard extrapolations of N-N collisions (also, we include corrections from incomplete thermalization of charm 
quarks in relaxation time approximation). Such an enhancement is needed if one aims at describing the NA50 data in 
terms of thermal production alone (assuming that directly produced charmonia have been completely dissolved). The 
NA50 suppression pattern towards more peripheral collisions then emerges as a canonical ensemble effect, requiring 
"anomalous" open-charm enhancement factors of around ^ 5 . Theoretically, such an increase in hard production 
in a heavy-ion environment is not easily justified. Also, the NA38 findings on the intermediate-mass dimuon excess 
allow for a maximal open-charm enhancement of a factor of ~3 for central collisions (less for peripheral)^ . Furthermore, 
imposing statistical production on small-size systems (or even p-A reactions) eventually contradicts nuclear absorption 
systematics. 

Finally we refer to Ref. where the J/'f abundance in an expanding QGP has been evaluated using rate 

equations with gluon photodissociation and its reverse reaction as microscopic input. The space-time evolution in 
ref. p8[ | occurs entirely within the QGP phase, and in-medium modifications to the J/^E" binding energy have not 
been applied. Consequently, the time evolution of J/^> formation in that approach is quite different from our results, 
leaning towards early production |Q, which also entails substantially larger results for the final Nj^^/Ncs ratio, 
reaching values of up to 20-30x10^'' for central Au-Au at y/sj^j^^— 200 GeV. This is a factor of 5-8 larger than our 
estimates, cf. Fig. ^ 



^ In our picture, this excess is attributed to thermal radiation |66l 
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VIII. CONCLUSIONS AND OUTLOOK 



In summary, we have developed a two-component model for charmonium production in heavy-ion colhsions within 
a comprehensive thermal evolution scenario. The two sources of charmonium production are (i) "direct" J/4''s, 
arising from primordial N-N collisions, subjected to nuclear, QGP and hadronic suppression, as well as (ii) statistical 
recombination of c and c quarks at the hadronization transition subjected to hadronic dissociation only. The QGP 
dissociation of contribution (i) has been evaluated using in-mcdium charmonium binding energies which led us to 
introduce "quasifree" destruction as the dominant suppression mechanism. Inelastic hadronic interactions have been 
estimated within earlier proposed effective lagrangian approaches, and turned out to give very moderate corrections 
(up to 10% for excited charmonium states). Contribution (ii) has been based on open-charm abundances as inferred 
from N-N collisions (without "anomalous" enhancement), and incomplete thermalization has been incorporated via 
a relaxation time approximation. Taken together, with an effective strong coupling constant as a single parameter, 
the measured centrality dependence of J/5' production at the SPS can be reasonably well described. A potential 
discrepancy has been identified in the ^E''/^' ratio, which we believe to reside in shortcomings of the calculations for 
hadronic ^' dissociation. The latter is difficult to assess in purely hadronic models and might well be underestimated. 

We have extrapolated our approach to higher energies. The pertinent excitation function for the Nj/x^/Ncc ratio 
in central collisions exhibits a non-trivial minimum structure around ^snn 40 GeV signalling the transition from 
(predominantly) direct to thermal production ||l^ . For the centrality dependence of this ratio at full RHIC energy, 
we predict a rather flat behavior for participant numbers beyond Npart — 150. The absolute number of produced J/^ 
mesons in central collisions at RHIC turns out to be close to what one expects from nuclear absorption alone, i.e., an 
80% QGP suppression is essentially regenerated by thermal production at hadronization. 

We emphasize again that the present analysis should also be considered as another step towards establishing 
a coherent picture of high-energy heavy-ion collisions as proceeding through a thermal evolution including QGP 
formation. Charmonium production has been linked with other observables such as hadro-chemistry dilepton 
production at low and intermediate mass |6^, 70 , etc. 

Concerning directions for future work, it will be necessary to reiterate the question of in-medium effects in hadronic 
dissociation cross sections especially for the lightly bound charmonium states. Here, quark- exchange models appear to 
be better suited than hadronic ones, allowing for modifications of effective confining potentials and constituent quark 
masses |Q - possibly related to effects of chiral restoration - on a microscopic level. Furthermore, the NA50 data on 
J/^E* transverse momentum distributions need to be addressed as another test of our approach. It is also mandatory 
to improve on our schematic treatment of the charm-quark thermalization. Here, valuable information can be expected 
from experiment in terms of both single- and di-lepton spectra from semileptonic open-charm decays 172, ^ . 
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